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ABSTRACT 

Three  methods  of  locating  events  with  regional  phases  were  tested  and 
compared  for  accuracy  of  hypocenter  determination.  A  modified  HYLO  method 
was  found  to  be  slightly  less  accurate  than  a  standard  location  algorithm 
using  a  Herrin  68  earth  model.  In  addition,  we  found  that  the  modified  HYLO 
method  requires  a  much  more  detailed  knowledge  of  the  earth's  crust  and  upper 
mantle  than  is  generally  available.  Two  methods  which  do  not  require  any 
previous  knowledge  of  the  crust  in  the  area  of  interest,  the  methods  of 
successive  determinations  and  simultaneous  inversions,  were  found  to  reduce 
location  errors  about  50%,  as  compared  to  the  standard  location  algorithm 
with  a  Herrin  68  earth  model. 
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1.  INTRODUCTION 


The  accuracy  of  location  of  events  located  with  regional  phases  (P  , 

P*,  P  ,  S  ,  S*,  S  )  is  affected  strongly  by  lateral  variations  in  crustal 
i  g  n  g 

and  upper  mantle  structures. 

Herrin  and  Taggart  (1962)  demonstrated  that  the  velocity  of  the  regional 

phase  P  varies  significantly  across  the  United  States.  Figure  1  shows  their 
n 

revised  map  of  apparent  P^  velocities  across  the  United  States.  Herrin  and 
Taggart  (1962)  introduced  a  method  known  as  HYLO,  which  applies  corrections 
to  the  travel  time  tables  of  phase  P^,  based  on  the  velocity  model  shown  in 
Figure  1.  They  were  able  to  improve  significantly  hypocenter  and  origin  time 
estimates  for  the  nuclear  explosion  GNOME  and  the  Hebgen  Lake  earthquake  of 
18  August  1959.  In  addition,  they  were  able  to  eliminate  the  bias  in 
travel  time  residuals  resulting  from  the  general  pattern  of  faster  P^ 
velocities  in  the  Eastern  United  States  (EUS),  as  opposed  to  the  Western 
United  States  (WUS).  Subsequently,  Herrin  and  Taggart  (1966)  reported 
similar  success  in  locating  the  nuclear  explosion  SALMON.  The  method,  however, 
has  several  drawbacks.  Only  the  phase  P^,  which,  in  many  regions  of  the  world, 
is  the  lowest  amplitude  phase  on  the  record,  may  be  used.  In  addition, 

HYLO  makes  no  correction  for  variation  in  crustal  thickness  between  source  and 
receiver,  which  as  Chang  and  Racine  (1979)  have  shown,  may  contribute  more 
to  travel  time  residuals  than  does  lateral  variation  in  the  phase  velocity. 

Chang  and  Racine  (1979)  relocated  twelve  nuclear  explosions,  utilizing 

data  for  the  phases  P  and  P  .  They  used  Julian's  (1974)  location  program  as 

a  g 

modified  by  McCowan  (1978)  to  accept  crustal  phases  P  ,  P*,  P  ,  S  ,  S*  and  S  . 

n  g  n’  g 

They  further  modified  the  program  to  allow  an  option  for  selecting  up  to 
fifty  local  crustal  models,  such  that  the  analyst  may  utilize  an  independent 
model  for  each  epicenter  to  station  path.  Figure  2  shows  the  boundaries  of 
the  regional  crustal  models  they  used,  first  published  by  Pakiser  and  Robinson 
(1966).  Figure  3  shows  the  boundaries  and  codes  for  fifteen  more  localized  models 
developed  by  Racine  (1979).  In  addition,  two  generalized  models,  Herrin  68  (1968) 
and  Jef f reys-Bullen  (1958)  were  also  used. 

Table  I  is  a  comparative  listing  of  the  general,  regional,  and  localized 
models.  Chang  and  Racine  (1979)  were  marginally  successful  in  improving 
the  location  estimates  of  the  twelve  explosions,  as  compared  to  a  standard 
location  algorithm  with  a  Herrin  68  model. 
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Figure 


TABLE  I 

General,  regional  and  local  crustal  models 


Model 

1st  Layer 

2nd  Layer 

Mantle 

Area 

Designator 

Thickness 

Velocity 

Thickness 

Velocity 

Velocity 

General: 

Jef frey-Bullen 

J-B 

15.00 

5.57 

18.00 

6.50 

7.80 

Herrin  68 

HE 

15.00 

6.00 

25.00 

6.75 

8.05 

Regional: 

Calif.  Coast 

ecus 

15.00 

6.20 

5.00 

7.00 

8.10 

Sierra  Nevada 

SNUS 

25.00 

6.20 

25.00 

7.00 

7.90 

Pac.  NW  Coast 

PCUS 

10.00 

6.20 

25.00 

7.00 

7.70 

Columbia  Plat. 

CAUS 

10.00 

6.20 

35.00 

7.00 

7.90 

Basin  &  Range 

BRUS 

20.00 

6.20 

10.00 

7.00 

7.90 

Colorado  Plat. 

CPUS 

25.00 

6.20 

15.00 

7.00 

7.80 

Rocky  Mtns. 

RMUS 

25.00 

6.20 

15.00 

7.00 

8.00 

Int.  Plains 

IPUS 

20.00 

6.20 

30.00 

7.00 

8.20 

Coastal  Plain 

CLUS 

20.00 

6.20 

15.00 

7.00 

8.10 

App.  Highland 
&  Sup.  Upland 

ASUS 

15.00 

6.20 

25.00 

7.00 

8.10 

Local: 

N.  Calif. 

NOCA 

12.00 

5.60 

18.00 

6.70 

8.00 

Coast  Calif. 

COCA 

10.00 

5.60 

10.00 

6.70 

8.00 

Sierra  Nevada 

SNCA 

15.00 

6.00 

20.00 

6.50 

7.60 

S.  Calif. 

SOCA 

20.00 

6.20 

10.00 

6.90 

7.80 

N.  Nevada 

NONV 

20.00 

6.00 

10.00 

6.70 

7.90 

Cent .  Nevada 

CENV 

20.00 

6.00 

10.00 

6. 60 

7.80 

SW  Nevdad 

SWNV 

27.00 

6.20 

9.00 

7.10 

7.80 

Lake  Mead 

Nev. 

LMNV 

15.00 

6.00 

15.00 

6.50 

7.90 

W.  Utah 

WEUT 

15.00 

5.90 

10.00 

6.40 

7.40 

E.  Utah 

EAUT 

27.00 

6.20 

13.00 

6.80 

7.80 

N.  Arizona 

NOAZ 

26.00 

6.00 

12.00 

6.80 

7.80 

Cent.  Arizona 

SOAZ 

15.00 

6.00 

7.00 

7.00 

7.80 

W.  Colorado 

WECO 

9.00 

6.10 

31.00 

6.60 

7.80 

W.  New  Mexico 

WENM 

19.00 

6.20 

21.00 

6.50 

7.90 
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Their  results  indicated  that  several  models  were  needed  for  each  event- 


station  path:  a  source  model  for  the  downgoing  ray,  a  path  model  to  give 
the  average  propagation  velocity,  and  a  receiver  model  for  the  upcoming  ray. 
The  original  experiment  allowed  only  one  model  per  path.  Their  results  also 


indicated  that  the  addition  of  P 

g 

location  accuracy. 


travel  times  did  not  significantly  improve 


This  year,  we  have  experimented  with  three  additional  methods  of  locating 

events  with  regional  phases.  The  first,  a  modified  HYLO  method,  incorporates 

an  option  for  selecting  crustal  models  at  the  source  and  receiver,  and  assumes 

a  constantly  dipping  Moho  between  the  two.  Two  other  methods  were  tested, 

which  require  no  previous  knowledge  of  crustal  structures,  and  which  may 

incorporate  phases  other  than  P  into  the  location  process.  The  method  of 

n 

successive  determinations  involves  alternately  calculating  a  least  squares 
fit  to  the  travel  times  of  the  arrivals  for  a  given  event,  and  then  relocating 
the  event.  Thus,  this  method  determines  an  independent  travel  time  table  for 
each  phase  of  each  individual  event.  The  method  of  simultaneous  inversions 
does  much  the  same,  except  that  the  new  travel  time  relationships  and 
locations  are  determined  simultaneously,  through  matrix  inversion,  as  opposed 
to  alternately,  as  in  successive  determinations. 


2.  THE  MODIFIED  HYLO  METHOD 


2.1  Theory  and  Method  of  Investigation 

The  original  HYLO  method,  of  Herrin  and  Taggart  (1962),  replaced  the 
travel  time  relationships  for  the  phase  from  the  standard  tables  with 
travel  times  calculated  by  determining  an  average  velocity  for  each  travel 
path,  using  a  velocity  which  varies  laterally  as  shown  in  Figure  1. 

Our  modified  HYLO  method  allowed  us  to  select  a  local  or  regional 
crustal  model  for  the  source  and  receiver  crust  from  among  those  shown 
in  Figures  2  and  3  and  Table  I.  A  constantly  dipping  Moho  was  assumed 
between  source  and  receiver.  Figure  1  was  divided  into  a  grid  of  one  degree 
squares,  defined  by  parallels  of  latitude  and  meridians  of  longitude.  Each 
block  of  the  grid  has  an  associated  P^  velocity.  A  computer  program  was 
developed  to  determine  a  theoretical  travel  time  for  a  ray  passing  down  the 
source  model,  across  the  grid  of  P^  velocities  along  the  constantly  dipping 
Moho,  and  up  the  receiver  model.  This  theoretical  time  was  subtracted  from 
the  expected  travel  time  taken  from  the  Herrin  68  table  for  the  same  distance, 
giving  us  a  station  correction  factor  to  the  Herrin  table,  which  is  valid 
for  the  given  path  and  crustal  models.  Since  we  were  using  data  from  nuclear 
explosions,  the  true  latitude  and  longitude  were  known  to  within  one 
hundredth  of  a  degree,  and  the  true  origin  times  were  known  to  within  one 
hundredth  of  a  second.  Therefore,  the  station  correction  factors  were  simply 
added  to  the  observed  P^  travel  times,  a  standard  location  program  was  run  on 
the  corrected  data,  utilizing  the  Herrin  68  travel  times,  and  absolute  errors 
were  obtained  by  comparing  the  computed  locations  with  the  true  values.  The 
results  were  compared  to  those  obtained  by  running  the  same  program  on 
uncorrected  data,  as  reported  by  Chang  and  Racine  (1979). 

2.2  Modified  HYLO  Method  -  Results  and  Discussion 

Table  II  presents  the  results  of  our  modified  HYLO  method.  The  average 
errors  in  location  and  origin  time  are  slightly  greater  than  those  obtained 
by  Chang  and  Racine  (1979)  for  the  same  events,  using  uncorrected  arrival  times 

for  P  only  and  both  P  and  P  ,  and  a  Herrin  68  earth  model. 

n  n  g’ 

The  poor  results  of  this  experiment  may  be  due  to  a  number  of  reasons. 

The  lateral  variation  in  P^  velocity  may  be  more  complicated  than  is  shown  in 
Figure  1.  Our  assumption  of  a  constantly  dipping  Moho  may  have  been  grossly 


Comparison  of  Location  Errors:  Herrin  68  Versus  Modified  HYLO 


inaccurate,  especially,  for  example,  when  a  phase  crosses  two  or  more 
distinct  geological  provinces,  as  in  a  path  from  the  Basin  and  Range  to 
the  Rocky  Mountains.  There  is  a  possibility  that  our  derived  local  models 
are  not  adequate.  We  find  a  great  deal  of  variation  among  authors  regarding 
local  crustal  models,  and  our  limiting  condition  that  the  crust  be  two  layers 
may  have  been  an  important  factor.  All  the  above  are  supported  by  Schilt  et 
al. (1979)  who  in  a  review  of  the  heterogeneity  of  the  crust,  with  data  from 
vibroseis  surveys,  carried  out  by  the  Consortium  for  Continental  Reflection 
Profiling,  conclude  the  following: 

1.  That  the  continental  crust  demonstrates  heterogeneities  on  a 
scale  of  a  few  kilometers  or  less,  as  well  as  patterns  on  a  continental 
scale. 

2.  That  the  Mohorovicic  discontinuity  has  the  characteristics  of  a 
complexly  layered  transition  zone,  as  opposed  to  a  continuous  interface 
separating  the  crust  and  mantle. 

Despite  all  those  possible  causes  of  error,  we  feel  that  our  revised  HYLO 
method  is  a  better  method  than  the  original  HYLO,  which  attempted  to  compensate 
only  for  the  portion  of  the  crust  above  sea  level.  We  find  that  the  HYLO  method, 
as  tested  in  the  previous  paper,  and  in  this  improved  version,  does  not  improve 
location  accuracy  if  only  a  few  stations  are  available  for  analysis;  in  Herrin's 
publications  the  number  of  stations  was  on  the  order  of  100  and  perhaps  local 
fluctuations  averaged  out  so  that  the  removal  of  overall  East-West  bias 
resulted  in  a  more  accurate  location. 
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3.  METHOD  OF  SUCCESSIVE  DETERMINATIONS 


I 


3.1  Review  of  the  P  and  P  Observations 

n  g 

Thus  far,  we  have  shown  that  the  modified  HYLO  method  was  not  effective 
in  eliminating  errors  in  location  due  to  the  earth's  heterogeneity.  We  have 
discussed  several  possible  explanations  of  why  it  did  not  work.  We  also 
discussed  the  fact  that  HYLO  is  difficult  to  use,  especially  for  smaller 
regional  events,  because  it  utilizes  only  the  phase  P  . 

The  reason  P  has  not  been  used  in  the  HYLO  method  is  that  it  has  been 
g 

commonly  understood  that  P  is  not  a  head  wave  and  that  P  is  not  linear  on 

8  8 

the  T-delta  curve.  Thus,  it  was  not  considered  adequate  for  use  in  locating 
events.  If  this  is  so,  it  may  also  be  possible  that  the  earth  is  so  hetero¬ 
geneous  that  P^  travel  times  fluctuate  widely  from  a  linear  relationship  so  that 
P^  may  not  be  suitable  for  locating  events.  We  have  investigated  this  point. 

Figure  4,  a  composite  travel  time  curve  for  eight  explosions  at  NTS, 

shows  that  the  travel  times  for  both  P  and  P  are  linear.  Thus,  we  conclude 

n  g 

that  both  P^  and  P  may  be  useful  for  location  purposes.  Figure  4  also  demonstrates 

that  it  is  possible  to  get  a  reasonably  good  average  velocity  for  phases  P^  and 

P  ,  thus  explaining  Chang  and  Racine's  (1979)  reasonably  good  results  in 
6 

locating  the  events  with  the  Herrin  68  earth  model.  Close  inspection  of  Figures 

4a  through  4n  and  Appendix  I,  however,  shows  that  while  the  travel  times  for  each 

individual  event  are  linear,  the  slopes  of  the  linear  fits  differ  for  sources 

distant  from  each  other,  and  in  no  case  are  the  velocities  exactly  equal  to  the 

Herrin  68  model.  Thus  we  conclude  that  we  may  improve  our  location  estimates, 

by  evaluating  P^  and  travel  time  relationships  for  each  event  and  its  unique 

station  distribution.  This  is  a  paradox,  of  course;  perhaps  the  P  and  P 

n  g 

velocities  are  characteristic  of  the  detailed  source  location  and  depth  and  not 
of  the  "broad-brush"  path  because  a  different  set  of  rays  (modes)  escape  the 
source  for  different  near-source  geologies. 

This  is  the  basis  for  our  methods  of  successive  determinations  and  simul¬ 
taneous  inversions.  We  are  not  trying  to  obtain  or  utilize  crustal  structures, 
which  we  have  demonstrated  may  vary  significantly  across  a  small  area.  We  are 

trying  to  determine  the  constant  P  and  P  velocities  that  best  fit  each 

n  g 

particular  event. 
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Figure  4a 


Observed  P  and  P  Travel  Tiroes  for  FAULTLESS 

n  8 
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Figure  4d  Observed  P  and  P  Travel  Times  for  ROCKVILLE  DAM 
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Figure  4f  Obse?-vcd  P  and  P  Travel  Times  for  KLICKITAT 
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Figure  4o  Observed  P_  and  P„  Travel  Times  for  NTS  COMPOSITE 


3.2  Theory  and  Method  of  Investigation 

The  method  of  successive  determinations  has  several  advantages  over  the 

modified  HYLO  method.  No  previous  knowledge  of  the  crust  is  necessary,  and 

it  is  then  applicable  anywhere  in  the  world.  The  method  is  applicable  to  any 

regional  phase,  whereas  HYLO  utilizes  only  the  phase  P^.  In  this  experiment, 

we  limited  ourselves  to  the  phases  P  and  P  ,  in  order  to  allow  ourselves 

n  g* 

to  compare  the  results  directly  with  Chang  and  Racine’s  (1979)  results. 

The  method  of  successive  determinations  utilizes  an  initial  epicenter 

location  which  may  be  determined  using  any  of  the  standard  travel  time 

tables.  The  location  used  in  this  manner  will,  of  course,  be  only 

approximately  correct,  and  the  size  of  the  absolute  error  of  location  will 

depend  on  the  difference  between  the  P  and  P  velocities  used  in  the  model 

n  g 

selected  and  the  true  velocities  for  these  phases  measured  along  the 

particular  source  to  receiver  paths  used  in  the  location.  In  this 

particular  study,  our  original  epicenter  determination  was  made  using  the 

Herrin  68  travel  time  tables  and  a  Herrin  68  crustal  model.  We  now  make 

the  critical  assumption  that  for  each  individual  event  the  P  and  P 

n  g 

velocities  have  constant  values  for  all  paths  and  hence  that  the  P  and  P 

n  g 

travel  time  "curves"  are  in  fact  straight  lines.  Evidence  substantiating  this 
assumption  has  been  presented  in  the  previous  section. 

We  shall  therefore  attempt  to  correct  the  approximate  locations  for  the 

errors  in  the  assumed  travel-time  relations,  but  in  practice  we  must,  of  course, 

do  so  without  any  reference  to  the  true  locations ,  which  were  known  a  p^LLOKi 

in  the  construction  of  Figures  4a-4o.  In  order  to  make  such  a  correction,  we 

construct  travel-time  curves  similar  to  those  in  Figure  4,  but  now  we  plot  the 

observed  arrival  times  versus  the  distance  from  the  approximate*  rather  than  the 

true,  epicenters.  A  linear  least-squares  fit  is  then  performed  on  these 

observations,  and  the  slopes  of  these  lines  yield  corrected  values  for  the 

P  and  P  velocities.  Next,  the  newly  determined  velocities  are  used  as  layer 
n  g 

velocities  in  a  crust  which  is  otherwise  appropriate  to  the  Herrin  travel 
times,  and  complete  travel  times  are  computed  by  raytracing.  This  process 
continues,  alternately  determining  phase  velocities  and  then  travel-times, 
location  and  origin  time,  until  the  solution  converges.  Thus,  the  method 
determines  a  new  set  of  travel  time  relationships  for  each  event. 

An  option  was  installed  to  allow  for  two  separate  sets  of  travel  time 
relationships  for  each  event.  This  option  may  be  useful,  for  example,  in 

-  -32- 


allowing  the  use  of  separate  travel  time  relationships  for  paths  crossing 
very  dissimilar  geological  provinces,  or  for  differentiating  between  stations 
on  opposite  sides  of  a  dipping  slab. 

3.3  Results  and  Discussion 


Table  III  presents  the  results  of  the  method  of  successive  determinations. 

Unlike  Chang  and  Racine’s  (1979)  earlier  findings  using  a  standard  location 

program  and  local  or  regional  crustal  models,  the  addition  of  P  arrival  times 

S 

to  the  data  base  of  the  method  of  successive  determinations  significantly 
improves  the  location  accuracy.  Thus,  it  may  be  that  the  tabulated  travel 


times  for  the  phase  P  are  not  applicable  to  wide  areas,  and  better  location 

S 

accuracy  is  obtained  by  determining  the  travel  time  relationships  for  P 

K 

for  each  individual  event.  This  may  be  due  to  the  fact  that  P  propagates 

S 

through  the  upper  layer  of  the  crust,  which  is  more  laterally  heterogeneous 


in  its  physical  properties  than  the  crust-mantle  interface  along  which  the 


P^  ray  propagates.  Or  it  may  be  possible  that  a  wider  variety  of  P^  rays  are 
emitted  as  a  function  of  source  depth  and  geology. 


Table  IV  compares  the  results  of  successive  determinations  versus  the 
location  and  origin  time  determined  using  a  standard  location  algorithm  and 
Herrin  68  earth  model.  Substantial  improvement  in  location  is  seen  using 
the  method  of  successive  determinations. 


Table  V  compares  the  results  of  successive  determinations  versus  the 

modified  HYLO  method.  Successive  determination  yields  a  more  accurate 

location,  requires  no  previous  knowledge  of  the  velocity  structure  of  the 

earth's  crust,  and  allows  the  use  of  phases  other  than  P  . 

n 

Table  VI  shows  the  results  of  the  double  successive  determination 
experiment.  This  experiment  invoked  the  option  of  separating  the  detecting 
stations  into  two  azimuthal  sectors  and  then  fitting  the  travel  times  for  each 
sector  separately.  As  is  shown  in  the  table,  dividing  the  arrival  times 
for  GNOME  into  eastward-going  and  westward-going  paths  does  reduce  the  loca¬ 
tion  error  somewhat.  However,  this  technique  fails  in  the  case  of  SHOAL,  as 
might  be  expected  since  all  directions  out  of  SHOAL  lie  within  the  Western  United 
States.  Note  that  the  velocity  for  the  second  group  of  stations  observing 
SHOAL  defaulted  to  the  Herrin  68  values,  when  the  least-squares  fit  to  the 
arrival  data  showed  too  much  variance.  Further,  for  SHOAL  the  first  seven 
iterations  appeared  to  diverge. 
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Figures  5a-n  show  the  actual  location  of  each  event  studied,  the  initial 
epicenter  determination  and  the  successive  locations  determined  by  our 
program  as  it  converges  on  a  solution.  The  open  circles  for  SHOAL  and  GNOME 
illustrate  the  successive  locations  determined  by  the  program  allowing  two 
separate  travel  time  relationships.  It  may  be  seen  that  most  solutions 
converged  after  about  seven  successive  approximations.  However,  the  event 
ROCKVILLE  DAM,  which  gave  a  good  solution  after  seven  successive  approxi¬ 
mations,  did  not  converge  to  a  final  solution. 
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it 

P  's 
n 

FAULTLESS 

5 

RULISON 

10 

PASSAIC 

4 

ROCKVILLE  DAM 

9 

DORMOUSE' 

3 

KLICKITAT 

13 

BANDICOOT 

5 

SHOAL 

12 

MERRIMAC 

8 

GASBUGGY 

12 

PILEDRIVER 

5 

ROANOKE 

4 

GNOME 

18 

TABLE  III 

Method  of  Successive  Determinations  -  Absolute  Errors 


Error  (km)  Error  (km)  Error  (km)  Error  (km) 
Pn  only  Pn  &  Pg  Pn  only  Pn  &  Pg 


it  P  'si  Depth  Free  Depth  Free  I  Depth  Rest.  iDepth  Rest. 


18.64 


(3.95) 


restr.  restr. 

(2.66)  restr. 

restr.  restr. 

diverges  (9.39) 

diverges  (3.11) 

N/A  restr. 

diverges  no  P 


restr.  no  P 

8 

10.51  restr. 

diverges  (1.83) 


18.64 

2.65 

2.66 

2.04 

diverges 

4.34 
N/A 

2.34 
5.17 

10.51 

10.41 


(4.41)  ignores  P  **  4.38 

g 

restr.  restr.  7.10 


2.14 
a  3.75 
b  5.00* 

10.70 
4.11 
9.50 
no  P 

g 

no  P 

g 

10.02 


Mean  Error  (km) : 


a  5.59 
b  5.71 


DIVERGES  =  Does  not  converge  to  solution 
RESTR  =  Restricted  to  0  depth 


Unlike  HYL0  -  can  use  P^  (bigger  than  P^,  especially  in  WUS) 
Unlike  Chang  &  Racine  -  P^  improves  location. 

*  Solution  alternates  between  two  epicenters 
**  Residuals  too  large  -  data  rejected 


TABLE  IV 


Accuracy  of  Location 

Successive  Determination  Versus  Herrin  68 
Depth  Free 


HERRIN  68 

SUCCESSIVE 

DETERMINATIONS 

LOC  E1 

OT  E2 

DEPTH3 

LOC  E 

OT  E 

DEPTH 

FAULTLESS 

15.04 

1.4 

0.0 

3.96 

1.0 

3.6 

RULISON 

1.58 

0.0 

0.0 

3.39 

2.6 

0.0 

PASSAIC 

4.08 

-0.2 

0.0 

2.14 

-0.6 

0.0 

5.18 

a  3.75 

1.9 

0.0 

ROCKVILLE  DAM 

1.5 

0.0 

b  5.00 

1.9 

0.0 

DORMOUSE' 

9.09 

0.2 

0.0 

9.87 

0.7 

4.3 

KLICKITAT 

6.53 

0.1 

0.0 

2.78 

0.7 

8.1 

BANDICOOT 

1.24 

0.6 

7.8 

4.58 

0.0 

0.0 

SHOAL* 

5.18 

0.0 

0.0 

2.34 

-0.9 

0.0 

MERRIMAC 

9.58 

0.3 

0.0 

5.16 

0.2 

15.9 

GASBUGGY 

11.43 

0.6 

0.0 

10.04 

2.0 

0.0 

PILEDRIVER 

12.98 

0.0 

0.0 

1.66 

-0.1 

8.9 

ROANOKE 

8.24 

1.3 

7.9 

4.20 

1.1 

28.1 

Average 

7.52 

a  4.49 
;  b  4.59 

*  Depth  restricted 

1.  LOC  E:  Error  of  location  vectors  in  km 

2.  OT  E:  Error  of  origin  time  estimates,  T  -  T  ,  in  sec 

est  true 

3.  when  depth  is  zero,  the  depth  is  restricted  to  0  km 

4.  Solution  alternates  between  two  epicenters. 
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TABLE  V 


Successive  Determination  Versus  HYLO 


HYLO 

LOC  E  OT  E 


SUCCESSIVE  DETERMINATION 
LOC  E  OT  E 


FAULTLESS 

12.32 

-1.8 

3.96 

1.0 

RULISON 

3.72 

1.0 

3.39 

2.6 

PASSAIC 

11.17 

0.8 

2.14 

-0.6 

ROCKVILLE  DAM 

5.59 

-0.4 

a  3.75 
b  5.00 

!-9  * 

1.9 

DORMOUSE' 

10.47 

-1.4 

9.87 

0.7 

KLICKITAT 

9.64 

0.6 

2.78 

0.7 

BANDICOOT 

7.09 

-1.8 

4.38 

0.0 

SHOAL 

11.26 

0 

2.34 

-0.9 

MERRIMAC 

10.89 

4.7 

5.16 

0.2 

GASBUGGY 

5.81 

0.4 

10.04 

2.0 

PILEDRIVER 

24.54 

0.9 

1.66 

-0.1 

ROANOKE 

3.68 

0.4 

4.20 

1.1 

Average 

9.68 

a  4.49 
b  4.59 

*  Solution  alternates  between  two  epicenters 


TABLE  VI 


Double  Successive  Determinations 


Group  //I 

GNOME 

Group  It  2 

Group  #1 

SHOAL 

Group  It 2 

LC-NM 

PO-TX 

EK-NV 

BMO 

ML-NM 

SS-TX 

WI-NV 

RT-NM 

SM-TX 

MV-CL 

CP-CL 

SV-AZ 

LP-TX 

CU-NV 

UBO 

SF-AZ 

HB-OK 

KN-UT 

BX-UT 

DR-CO 

AM-OK 

HL-ID 

TFO 

FS-AZ 

TO-OK 

DR-CO 

WM-AZ 

SJ-TX 

KN-UT 

PM-WY 

7.887' 

8.244 

8.108 

*** 

LOC 

Error 

SUCCESS 

(km)  OT  Error  (sec) 

H  (km) 

DOUBLE  SUCCESS 

LOC  Error  (km)  OT  Error  (sec) 

H  (km) 

GNOME 

7.07 

+0.3 

0.0* 

5.73  +1.5 

4.965 

SHOAL 

2.34 

-0.9 

0.0* 

11.04 (diverges)  +0.9 

0.0* 

*  RESTRICTED 


***  DEFAULTED  to  Herrin  68  Model 


FAULTLESS 

o 


Figure  5a  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 


Figure  5b  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 


116  080  116.060  116.040 

Figure  5c  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 


106.460  106.440  106.420  106.400 

Figure  5d  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 


Figure  5e  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations 


KLICKITAT 


Figure  5f  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations 


Figure  5k  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 


Figure  51  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 


Figure  5m  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 
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Figure  5n  Convergence  of  calculated  epicenters  toward  true 
epicenters  -  Method  of  successive  determinations. 
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4.  METHOD  OF  SIMULTANEOUS  INVERSIONS 


4.1  Theory 

Thus  far  in  this  report  two  approaches  have  been  taken  to  the  problem 
of  overcoming  inaccuracies  introduced  into  the  location  of  regional  events 
by  an  inadequate  knowledge  of  the  and  travel-time  relations.  The 
first  approach,  HYLO,  involved  the  construction  of  improved  travel-time 
curves  using  velocities  measured  by  crustal  refraction  methods.  The 
second  approach,  the  method  of  successive  determinations,  assumes  no  such 
a  pnj.o>vL  information  (except  as  an  initial  guess),  but  instead  determines 
corrections  to  the  assumed  travel-time  curves  alternately  with  corrections  to 
the  assumed  event  locations  in  a  successive  approximations  scheme.  We  shall 
now  describe  a  third  approach,  the  method  of  simultaneous  inversions,  whereby 
the  corrections  to  the  travel-time  curves  and  to  the  event  locations  are 
computed  at  the  same  time  rather  than  alternately.  We  shall  present  a 
detailed  comparison  of  these  latter  two  methods  in  order  to  determine  which 
approach  ought  to  be  taken  in  order  to  improve  the  location  of  events  at 
regional  distances  in  the  absence  of  CL  psiiotii  knowledge  of  the  true  travel¬ 
time  relations. 


The  standard  technique  for  locating  seismic  events,  at  both  regional 

and  teleseismic  distances,  involves  assuming  a  trial  value  of  the  hypocenter 

and  origin  time,  calculating  predicted  arrival  times  at  a  set  of  stations  for 

signals  originating  at  the  trial  hypocenter,  expanding  the  travel  time 

residual  at  the  i^  station,  6  t .  =  t,  ,  ,  -  t.  ,  n  _  .»  in  the  first- 

i  i, observed  i,  calculated 

order  Taylor  series 


5ti  •  #  >1 dT  +  <£>1 d*  +  <F>i dy  +  (3t>i  dh- 


(1) 


forming  a  system  of  such  equations  for  all  residuals  6  t^,  solving  for  the 
location  corrections  dT,  dx,  dy,  and  dh  by  means  of  matrix  inversion  and  then 
adding  these  corrections  to  the  assumed  values  of  the  event's  origin  time,  east, 
north,  and  depth  coordinates.  In  order  to  evaluate  the  derivatives  in  this 
equation,  one  may  set 

(2) 


at 

—  1 

3T 

B  X 

at 

at 

3A 

3x 

"  a  a 

3x 

at 

at 

3A 

3y 

”  3A 

9y 

^  ,  at 

sin  co  *  JK 


r  3t 

COS  C0  *  3A 


(3) 

(4) 


-46- 


■Afchriyfti  n  a.-fifrr 


where  is  the  azimuth  from  the  epicenter  to  the  station.  If  the  travel- 

8 1  8 1 

time  relations  are  accurately  known,  then  the  derivatives  —  and  —  may  be 

calcuated  at  the  distance  and  the  matrix  inversion  thereby  made  tractable. 

Evidence  has  been  presented  in  Figure  4  and  Table  III  that  the  travel-time  curves 

for  P  and  P  are  linear,  and  we  shall  therefore  assume  that  the  travel  times 
n  g 

may  be  defined  analytically  as 


tp  =  T  +  ap  (h)  +  bp  -A  (5) 

n  n  n 

tp  =  T  +  ap  (h)  +  bp  -A.  (6) 

g  g  g 


In  these  expressions  we  have  made  it  explicit  that  the  coefficients  ap^  and  ap^ 
depend  on  the  depth  of  the  event,  since  they  represent  the  delay  between  the 
event  origin  time  and  the  first  motion  at  the  surface.  The  travel-time 
residuals  at  regional  distances  may  thus  be  expanded  as 


©t  )±  =  Otp  /3T).  dT  +  (3t  /3x).  dx  +  (3tp  /3y)  dy 

n  n  n  n 

+  (3tp  /3ap  )i  dap  +  Otp  /9bp  >.  dbp 
n  n  n  n  n  n 

(5  t  )  =  (3tp  / 3T)  dT  +  (3tp  /3x).  dx  +  (3tp  /3y).  dy 

8  g  8  8  1 

+  Otp  /3ap  )±  dap  +  (3tp  /3bp  )±  dbp  . 
g  g  g  g  g  g 


(7) 


(8) 


These  expressions  are  not  explicitly  a  function  of  depth,  since  the  depth 

dependence  is  contained  within  the  coefficients  a  and  a  .  Evaluation  of 

*  n  *  g 

the  derivatives  in  these  expressions  is  straightforward: 


9 1  /3T  =  3t_  /3T  =  3tD  OaP„  =  3tP  /3aP  =  1 

pn  pg  pn  1  g  pg 

3tp  /3x  =  -  sin  •  (3tp  /3 A)  =  -  sin  •  bp 

n  n  n 

9tp  /3x  =  -  sin  •  bp 

g  g 

3tp  /3y  =  -  cos  •  (3 tp  /3A)  *  -  cos  •  bp 


(9) 

(10) 
(11) 
(12) 
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9tp  / 3y  =  -  cos  •  bp  (13) 

g  g 

Otp  /3bp  )  -  A±  (14) 

rn  n  i 

(3tp  /3bp  )±  =  A1.  (15) 

g  g 

The  addition  of  the  travel-time  residuals  5  t  and  5  p  to  the  conventionally 

n  g 

calculated  residuals  5  t„  ,  .  enables  the  event  location  to  be  expressed 

P  teleseismic 

in  terms  of  an  expanded  set  of  coordinates  (x,y,h,T,ap  ,bpn,apg,bpg) .  For  an 
over-determined  system  of  more  than  eight  observations,  the  matrix  inversion 
may  be  carried  out  as  before  and  a  solution,  in  the  least-squares  sense,  may 
be  found  for  the  newly  defined  eight-dimensional  hypocenter.  The  first  four 
coordinates  give  the  event  location,  and  the  last  four  determine  the  slopes  and 
intercepts  of  the  and  Pg  travel-time  curves.  If  there  are  no  observations  of 
teleseismic  P,  then  the  system  of  equations  does  not  involve  h  explicitly  and  so 
the  depth  cannot  be  determined,  but  the  other  seven  hypocentral  coordinates  still 
can  be.  Similarly  the  absence  of  P^  or  P  arrival  time  data  presents  no  dif¬ 
ficulties  in  the  computation  since  the  matrix  manipulation  is  simply  carried  out 

in  a  six-dimensional  subspace  by  restraining  the  unmeasured  coefficients  aP 

rn 

and  bp  or  ap  and  bP  to  have  some  arbitrary  fixed  value. 

*n  rg  rg 

In  the  event  that  no  teleseismic  P  arrival  times  are  available  and  that  all 

of  the  regional  data  consist  of  measurements  made  upon  a  single  phase  only  (i.e., 

all  P  or  all  P  ) ,  then  the  origin  time  as  well  as  the  depth  of  the  event 
n  8 

becomes  indeterminate.  That  this  is  so  can  be  seen  by  examining  equations  (5) 

and  (6) ,  wherein  it  is  shown  that  in  the  absence  of  P  measurements  one  may 

calculate  only  the  sum  of  T  and  aP  rather  than  the  two  individual  coefficients, 

rn 

with  a  similar  indeterminacy  occurring  in  the  case  of  no  P^  data. 

A  consideration  of  the  physical  meaning  of  the  coefficients  ap  ,  bp  ,  ap  , 
and  bp  permits  the  depth  to  be  estimated  when  it  cannot  be  directly  calculated 
on  accSunt  of  the  absence  of  teleseismic  data.  Consider  a  simple  earth  model 
consisting  of  an  upper  crustal  layer  of  thickness  h^  and  P-wave  velocity,  V^, 
a  lower  crustal  layer  of  thickness  and  velocity  V^,  and  a  mantle  with 
velocity  V^.  The  ray  diagram  in  Figure  6  shows  that 

tp  =  h/2VL  +  A/V  ,  (16) 

g 


so  we  may  set 


pr 


■w 


ap  =  h/2V 

(17) 

g 

bp  =i/vr 

(18) 

8 


In  a  like  manner  we  find  that  for  an  event  occurring  in  the  upper  crustal 
layer  (cf.  Figure  7a) 


ap  =  2hp/Vp  cos  0.^  +  2h2/V2  cos  02  -  (2^  tan  0p  +  2h2  tan  Q2^V3 
n 

+  (tan  0^/V^  “  1/Vp  COS  ®1^  *  k  (19) 


and 


for  an  event  occurring  in  the  lower  layer  (cf.  Figure  7b) 

ap  =  (hp  +  2h2)/V2  cos  0 2  +  h^/V^  cos  0p  -  [  (hp  +  2h2)  tan  ©2  +  hp 
n  +  (tan  92/V3  "  l/v2  cos  02)  *  h 


tan  011/V3 
(20) 


and  in  both  cases 


bp  =  1/V3  (21) 

n 

where  the  angles  in  these  expressions  are  given  by  Snell's  law  as 


sin  0p  =  Vp/V3 

(22) 

sin  ©2  *  V2^V3’ 

(23) 

We  shall  abbreviate  the  preceding  expressions  for  ap  and  ap 

p 

and  lower  crustal  layers)  as  6 

(for  both  the  upper 

aP  =  ki  *  h 

(24) 

g 

aP  =  k2  *  h  +  k3 

(25) 

n 

and  form  the  difference 

aP  “  aP  =  (kl  "  k2^  *  h  "  k3' 

(26) 

g  8 
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Substituting  into  the  above  formula  the  values  which  are  found  for  ap  and 

S 

Up  by  the  method  of  simultaneous  inversions,  we  may  now  calculate  the  depth 
innterms  of  the  model  parameters  k^,  k£,  and  k^.  We  emphasize  that  such  a 
calculation  is  only  approximate,  since  the  impetus  for  using  the  simultaneous 
inversion  approach  is  the  assumption  that  the  a  p^viofu.  earth  model  is 
inaccurate.  Equations  (24)  and  (25)  cannot  be  used  to  determine  the  depth  if 

only  P  or  P  data  are  used,  since  in  these  cases  only  T  +  a_  or  T  +  a„  can 
be  determined.  If,  however,  the  location  is  constrained  to  tRe  surf ace ,pthen 
these  expressions  show  that  the  origin  time  may  be  approximated  by  means  of 
relations 


ap  (measured)  «k  *h  +  T  =  T  (27) 

g 

ap  (measured)  =  (k^  -  k^)  •  h  -  k^  +  T  =  -k^  +  T.  (28) 

n 

Table  VII  presents  a  summary  of  which  variables  can  be  determined  when 
different  types  of  data  are  used.  The  indeterminate  variables  represent  a 
difference  between  the  method  of  simultaneous  inversions  and  the  method  of 
successive  determinations.  Even  though  both  methods  are  based  upon  the 
determination  of  the  slopes  and  intercepts  of  the  P^  and  P^  travel-time  curves 
and  hence  they  both  should  be  indeterminate  under  the  same  circumstances,  differ¬ 
ences  in  the  operational  approaches  which  were  taken  to  the  implementation  of 
the  two  techniques  create  discrepancies  in  their  applicability.  The  simultaneous 
inversions  technique  is  based  directly  on  equations  (5)  and  (6),  which  cause  the 
indeterminations  shown  in  Table  VII;  the  implementation  of  the  successive 
determinations  method,  however,  uses  a  hybrid  approach  which  calculates  the  slopes 
and  intercepts  in  equation  (5)  and  (6)  but  then  ignores  the  intercepts  and  substitutes 

the  resulting  velocities  in  a  Herrin  structure  which  may  be  used  to  calculate 
dT/dA  and  dT/dh  for  regional  phases,  equations  (17)  -  (23).  During  each 
iterative  calculation  of  the  hypoc'nter,  the  travel-time  curves  are  thus  taken 
to  be  pre-determined ,  so  the  indeterminations  shown  in  Table  VII  do  not  occur. 

In  the  comparison  of  results  of  the  two  methods  which  will  be  presented  in  the 
next  section,  there  are  therefore  several  instances  of  the  calculation  of  some 
variable  by  the  successive  determinations  method  even  though  no  corresponding 
value  can  be  calculated  by  the  simultaneous  inversions  method. 
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TABLE  VII 

Variables  Which  May  Be  Determined  Using  the  Method  of  Simultaneous  Inversions 


Determinate  Variables 


P  P  P 

n’  %'  teleseismic 


x,y,h,T,a  ,  bp  ,  ap  ,  b 

i  _  n  *  r 

n  g  g 


p  p 

n’  teleseismic 


x,y,h,T,ap  ,  bp 
n  n 


P  P 

g’  teleseismic 


x,y,h,T,ap  ,  bp 
g  g 


teleseismic 


inapplicable 


P  ,  P 
n  g 


x»y>T»3p  »bp  *^p 
n  n  g  g 
h  (by  approximation) 


x,y,bp  ,  (T  +  ap  ) 
n  n 


T  (by  approximation  if  depth  is  restrained) 


x,y,bp  ,  (T  +  ap  ) 
g  g 

T  (by  approximation  if  depth  is  restrained) 


cor  events  witn  epicenters  close  to  each  other  and  at  the  same  depth,  tin 
coefficients  ,  bp  ,  ap  ,  bp  should  be  nearly  constant  from  event  to  event. 
The  accuracy  of^the  §eterminat?on  of  these  coefficients,  and  hence  of  the 
event  location,  may  therefore  be  improved  by  inverting  the  arrival-time 
measurements  from  several  nearby  events  simultaneously,  a  process  known  as 
joint  epicenter  determination  (JED).  Applying  JED  to  data  from  N  events  has 
the  effect  of  transforming  N  systems  of  equations  in  8N  unknowns  to  a  single 
system  (with  the  same  total  number  of  equations)  in  4N  +  4  unknowns.  Extension 
of  conventional  location  to  encompass  JED  is  a  straightforward  technique 
described  by  Douglas  (1967)  and  by  Ahner,  Blandford  and  Shumway  (1971).  We 
shall  examine  whether  this  same  extension  can  be  performed  for  the  method  of 
simultaneous  inversions. 


4.2  Results 


In  order  to  Implement  the  method,  the  standard  location  program  LOCATION 
was  rewritten  to  solve  for  eight,  rather  than  four,  unknowns.  The  revised 
program  was  tested  using  the  same  data  set  as  was  used  for  testing  the  method 
of  successive  determinations.  (There  were  a  few  small  changes  from  the  data 
used  in  the  previous  section  of  this  report.)  Table  VIII  presents  the  results 
of  eight  trials  which  should  be  examined  in  comparing  the  two  methods.  The 
first  trial  represents  the  "correct"  results,  namely  the  travel-time  curves 
which  are  constructed  as  least-squares  fits  using  the  known  hypocenter  and 
origin  time,  as  shown  in  Figure  4.  Next  are  presented  the  results  of 
applying  the  method  of  successive  determinations  for  data  only  and  for  combined 
and  P^  data.  These  two  trials  are  then  repeated  with  the  depth  restrained  to 
the  surface.  The  simultaneous  inversions  results,  which  follow,  ought  to  be 
compared  with  those  obtained  by  the  successive  determinations  techniques  run  in 
both  the  depth-free  and  depth-restrained  modes.  The  reason  for  this  ambiguity 
is  that  even  though  the  method  of  simultaneous  inversions  was  run  in  the  depth- 
restrained  mode  (since  no  teleseismic  data  were  used) ,  the  depth  does  neverthe¬ 
less  vary  indirectly  through  the  coefficients  ap  and  ap  ,  hence  the  results  are 
perhaps  more  nearly  analogous  with  those  obtaineS  for  th§  depth-free  mode  of  the 
other  technique. 


It  was  not  possible  in  every  case  to  run  the  successive  determinations  program 
successfully  in  the  depth-free  mode.  The  depth  would  frequently  turn  out  to  be 
negative,  in  which  case  it  would  be  restrained  to  the  surface;  no  results  are 
given  for  these  cases,  since  they  are  identical  to  those  which  are  produced  by 
the  running  of  the  program  in  the  depth-restrained  mode  directly.  Some  other 
cases  resulted  in  negative  depths  which  were  restrained  to  the  surface  at  the 
end  of  the  first  iteration,  but  they  yielded  positive  depths  when  the  improved 
travel-time  curves  were  applied  to  them  in  the  second  iteration.  These  once- 
restrained  cases  are  denoted  by  parentheses  placed  around  the  value  of  the 
absolute  error. 


There  are  many  gaps  in  the  values  given  in  Table  VIII  for  the  method  of 
simultaneous  inversions  which  reflect  the  indeterminations  listed  in  Table  VII. 

As  has  been  discussed  previously,  the  inde termination  of  the  depth  can  be  removed 
by  approximation  if  one  assumes  an  earth  model  such  as  is  shown  in  Figures  6  and  7. 
We  have  assumed  the  following  parameters  for  such  a  model: 
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Figure  7a  Raypath  for  P  . 


Theoretical  Travel  Time  for  in  o  Simple  Earth  Model 
(2)  Hypocenter  in  Lower  Crustal  Layer 


X(S  (h|+h2-h)tan02 

S{ s  Xi/sin02 

X3  =  h2  tan  62 

S3  s  X3/sin^ 

X4  ~  h,tan0i 

S4  -  X4/sin0i 

S2  s  A-X,~  X3- X4 

tp„  *  Sj  /  V2  +  S2  /  V3 + S3  /  V2  +  S4  /  \^ 

=  (h|  ■►2h2 — h)/\^2  cos^* 

[A-  (h,+2h2  -  h)  ton  $>-  h,  tan  0j/ V- 

♦hi/\^C08®| 

Figure  7b 

Raypath  for  P  . 

TABLE  VIII 


Comparison  of  Method  of  Simultaneous  Inversions 
with  Method  of  Successive  Determinations 


Trial 

Trial 

Trial 

Trial 

Trial 

Trial 

Trial 


Trial  8 


least-squares  fit  of  travel  times  from  the  known  epicenter 
successive  determinations,  Pn  data  only;  depth  free 
successive  determinations;  Pn  and  Pg  data;  depth  free 
successive  determinations;  P 
successive  determinations;  Pn  and  P 
simultaneous  inversions;  Pr 
simultaneous  inversion;  Pn  and 


n  data  only;  depth  restrained 

g  data;  depth  restrained 
data  only 


Pg  data 


simultaneous  inversion;  Pn  and  Pg  data;  joint  epicenter  determination 

Vp 


Depth 

Absolute 

L 

n 

Degrees  of 

Event  Trial  // 

Lat(N) 

Long(W)  (km) 

0.  T.  Error  (km)(km/sec)(km/sec 

)  X  Freedom 

FAULTLESS 

1 

38.634 

11.216 

18:15:00.1 

5.882 

7.719 

0.64,  1.61 

3,3 

2 

.524 

.378  7.8 

:02.0 

18.64 

8.089 

1.72 

1 

3 

.625 

.260  3,6 

:01. 1 

(3.95) 

5.967 

7.922 

1.26 

6 

4 

.524 

.378 

:01. 7 

18.64 

8.088 

1.72 

2 

5 

.628 

.262 

:00.9 

4.05 

5.967 

7.922 

1.30 

7 

6 

.433 

.217 

22.35 

7.640 

1.26 

4 

7 

.468 

.212  -4.2 

:  59. 8 

18.46 

5.834 

7.648 

1.26 

4 

8 

.471 

.247 

:01 . 6 

18.31 

RULISON 

1 

39.406 

107.949 

21:00:00.1 

6.111 

8.192 

2.67,  14.06 

5,8 

2 

[depth 

restrained ] 

3 

[depth 

restrained ] 

4 

.421 

.937 

:  01 . 2 

2.65 

8.161 

1.91 

7 

5 

.398 

.987 

:  02 . 7 

3.38 

6.169 

8.292 

8.62 

14 

6 

.420 

.939 

1.78 

8.177 

13.63 

6 

7 

.412 

.915  -7.8 

:00. 3 

0.69 

6.111 

8.184 

16.44 

11 

8 

[too  far] 

PASSAIC 

1 

37.120 

116.040 

18:00:00.2 

6.536 

7.812 

0.04,  0.16 

1,2 

2 

.120 

.070  0.9 

:  59. 7 

(2.66) 

7.910 

0.00 

0 

3 

[depth 

restrained ] 

4 

.120 

.070 

:  59. 6 

2.66 

7.909 

0.00 

l 

5 

.122 

.064 

:  59. 6 

2.14 

restr . 

7.881 

1.04 

4 

6 

.120 

.070 

2.66 

7.907 

0.00 

0 

7 

.120 

.071  17.6 

:06.0 

2.75 

6.630 

7.909 

0.00 

1 

8 

.133 

.039 

:00. 7 

1.45 

ROCKVILLE 

DAM 

1 

39.360 

106.460 

16:21:33.6 

5.900 

8.213 

6.27,  6.80 

6,7 

2 

[depth 

restrained ] 

3 

.361 

.436  12.9 

:  37. 5 

(2.07) 

6.129 

8.139 

16.15 

13 

4 

[depth 

restrained ] 

5a 

.345 

.421 

:  35. 5 

1.90 

5.813 

8.127 

19.25 

14 

5b 

.316 

.426 

:  35 . 5 

1.90 

5.884 

8.136 

7.99 

14 

6 

.372 

.440 

2.59 

7.800 

2.81 

5 

7 

.336 

.419  -31.9 

:  27. 9 

4.42 

5.808 

8.102 

8.86 

1 1 

8 

[too  far] 
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TABLE  VIII  (con't) 


Depth 

Absolute 

VP„ 

VP 

1  n 

Degrees  of 

Event  Trial  It 

Lat(N) 

Long(W)  (km) 

0.  T.  Error  (km)(km/sec)(km/sec] 

1  x2 

Freedom 

DORMOUSE' 

1 

37.040 

116.020 

18:00:00.1 

6.234 

7.829 

0.39,  0.04 

2,1 

2 

[fails 

to  converge] 

3 

.116 

.974  4.3 

:00. 8 

(9.39) 

6.134 

7.969 

0.18 

3 

4 

[ fails 

to  converge] 

5 

.126 

.966 

:00. 5 

10.70 

6.119 

7.985 

0.19 

4 

6 

[insufficient  data] 

7 

.111 

.993  -6.7 

:00.4 

8.25 

6.159 

7.995 

0.24 

1 

8 

.077 

.994 

:00.8 

4.72 

KLICKITAT 

1 

37.151 

116.040 

15:30:00.1 

6.091 

7.762 

1.26,  14 

.88 

10,11 

2 

[fails 

to  converge] 

3 

.158 

.074  8.1 

:00. 8 

(3.11) 

6.108 

7.935 

6.22 

21 

4 

.152 

.089 

:  59.6 

4.34 

7.916 

3.31 

10 

5 

.158 

.089 

:00.2 

4.41 

6.112 

7.891 

8.01 

22 

6 

.146 

.031 

0.97 

7.900 

15.24 

9 

7 

.148 

.034  2.5 

:00.9 

0.63 

6.088 

7.759 

16.55 

19 

8 

.149 

.027 

:  00. 2 

1.18 

BANDICOOT 

1 

37.040 

116.020 

18:00:00.0 

5.763 

8.101 

0.17,  0.38 

2,3 

2 

.077 

.012  0.0 

:00. 3 

4.18 

8.205 

0.25 

1 

3 

[depth  restrained] 

4 

.077 

.012 

: 00.  3 

4.18 

8.205 

0.25 

2 

5 

.079 

.014 

:00. 5 

4.37 

6.115 

8.217 

2.40 

6 

6 

.140 

.013 

11.14 

7.900 

0.03 

1 

7 

.118 

.021  -18.9 

:  58. 9 

8.67 

5.983 

8.453 

0.32 

3 

8 

.995 

.012 

:00.1 

5.05 

SHOAL 

1 

39.200 

118.380 

17:00:00.1 

no  data 

7.831 

12.02 

10 

2 

[fails 

to  converge] 

3 

lno  Po 

data] 

4 

.194 

.406 

:  59. 2 

2.34 

7.893 

2.92 

3 

5 

[no  P„ 

data] 

6 

.187 

.394 

1.88 

7.835 

12.02 

8 

7 

[no  PK 

data] 

8 

[too  far] 

MERRIMAC 

1 

37.055 

116.033 

16:00:00.1 

no  data 

7.689 

2.62 

6 

2 

[depth  restrained] 

3 

[no  P„ 

data] 

4 

.01 J 

.058 

:  58. 6 

5.17 

7.731 

1.10 

5 

5 

[no  P 

data] 

6 

.01§ 

.058 

5.17 

7.731 

1.44 

5 

7 

[no  Pg 

data] 

8 

.01Z 

.059 

:  58.3 

5.31 

CAS BUGGY 

1 

36.678 

107.308 

19:30:00.1 

6.166 

8.260 

4.90,  13 

.88 

10,10 

2 

.635 

.203  -0.1 

:02. 3 

10.51 

8.203 

2.71 

8 

3 

[depth  restrained] 

4 

.635 

.203 

:02. 3 

10.51 

8.203 

2.71 

9 

5 

.  644 

.204 

:02.0 

10.02 

5.940 

8.198 

17.73 

21 

6 

.634 

.203 

10.56 

8.208 

2.07 

8 

7 

.633 

.210  -5.5 

:01.9 

10.07 

6.142 

8.214 

5.78 

18 

8  [too  far] 
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TABLE  VIII  (con’t) 


/  Lat(N) 

Depth 

Absolute 

VP 

e 

VP 

n 

„  Degrees 

Event  Trial 

Long(W)  (km) 

0.  T .  Error  (km)(km/sec)(km/sec) 

X 

Freedom 

PILEDRIVER  1 

37.227 

116.056 

15:30:00.1 

6.128 

7.729  0. 

29,  1.11 

3,3 

2 

[fails 

to  converge] 

3 

.223 

.076  8.9 

:00. 1 

(1.83) 

6.095 

7.7AO 

0.67 

6 

A 

.153 

.128 

:59.0 

10. A1 

7.823 

0.22 

2 

5 

.217 

.082 

:59.3 

2.56 

6.095 

7.7AO 

1.02 

7 

6 

.165 

.118 

8.81 

7.789 

0.A6 

1 

7 

.190 

.100  0.9 

:  00. 6 

5.67 

6.158 

7.76A 

1.01 

A 

8 

.179 

.090 

:00.3 

6.13 

ROANOKE  1 

37.123 

116.051 

15:00:00.2 

5.788 

7.A9A  0. 

AA,  0.29 

2,2 

2 

.160 

.033  2.5 

:59.2 

(A.A1) 

7.576 

0.02 

0 

3 

[drops 

P„  on  account  o 

f  depth] 

A 

.160 

•  03A 

:58.9 

A. 38 

7.576 

0.02 

0 

5 

.156 

.027 

:59.0 

A.2A 

5.893 

7.557 

3.A8 

5 

6 

.156 

.033 

A.  00 

7.539 

0.00 

0 

7 

.153 

.  02A  0.9 

:59.9 

A. 11 

5.81A 

7.529 

0.18 

2 

8 

.186 

.020 

:01 .0 

7.52 

SALMON  1 

31.1A2 

89.570 

16:00:00.0 

no  data 

8.263 

1.61 

3 

2 

[ fails 

to  converge] 

3 

[no  P 

data] 

A 

.  A26 

.671 

:01.1 

33.00 

8.138 

0.31 

2 

5 

[no  P 

data] 

6 

30.211 

.  2A6 

107.75 

8.290 

0.18 

1 

7 

[no  P_ 

data] 

8 

[too  far] 

GNOME  1 

32.263 

103.865 

19:00:00.0 

1 

data  pt 

7.920 

20.80 

16 

2 

[depth  restrained] 

3 

[depth 

restrained ] 

A 

.239 

.795 

:00. 3 

7.10 

7.983 

5.81 

15 

5 

.238 

.796 

:00. 3 

7.06 

restr . 

7.988 

5.90 

16 

6 

.238 

.795 

1 .15 

7.993 

9.02 

1A 

7 

.238 

.795  -0. A 

:  00. 5 

7.15 

6.095 

7.993 

9.02 

13 

8 

[  too 

far] 
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=  20.0  km 
=  10.0  km 
=  6.20  km/sec 
=  7.00  km/sec 
=  7.90  km/sec. 

Some  idea  of  the  validity  of  this  model  and  hence,  of  the  depth  approximation, 
may  be  gained  by  comparing  with  the  value  listed  for  the  velocity  of  P  , 

Vp  =  1/bp  ,  and  by  comparing  V  with  Vp  =  1/bp  . 

rg  &  j  rn  n 

The  absolute  errors  shown  in  Table  VIII  are  summarized  in  Table  IX.  It  is 
apparent  that  the  addition  to  the  data  base  of  P  arrival  times  frequently 
improves  the  location  relative  to  that  which  would  be  determined  using  P^  data 
alone,  and  it  seldom  worsens  it.  This  improvement  holds  for  both  the  method 
of  successive  determinations  and  the  method  of  simultaneous  inversions,  and  it 
represents  an  important  departure  from  the  results  obtained  by  conventional 
location  techniques  (Herrin  and  Taggart,  1962;  McCowan  and  Needham,  1978; 

Chang  and  Racine,  1979).  A  comparison  of  the  results  obtained  by  the  method 
of  successive  determinations  and  by  the  method  of  simultaneous  inversions 
shows  that  even  though  one  or  the  other  of  the  two  methods  yields  a  better 
location  for  particular  events,  on  the  whole  the  results  are  about  the  same. 

A  notable  exception  is  SALMON,  for  which  the  method  of  successive  approxima¬ 
tions  yielded  a  poor  result  and  the  method  of  simultaneous  inversions  failed. 
This  failure  is  due  to  the  fact  that  the  least-squares  solution  had  only  one 
degree  of  freedom  (cf.  Table  VIII).  The  fact  that  there  are  discrepancies 
between  the  two  methods  suggests  that  in  practice  perhaps  both  techniques 
should  be  used  in  the  location  of  unknown  events;  subjective  judgment  would 
then  be  necessary  to  select  one  of  the  two  resulting  epicenters  as  "the" 
location. 

Tables  VIII  and  IX  also  show  the  results  of  applying  JED  to  those  eight 

NTS  events  for  which  the  composite  travel-time  curves  are  shown  in  Figure  4o. 

As  one  might  expect  for  a  technique  such  as  this  which  effectively  performs 

an  average  over  events,  in  some  cases  the  location  was  improved  but  in  other 

cases  it  was  worsened.  The  net  effect  seems  to  be  small.  No  depths  are 

approximated  for  the  trial  using  JED  since  this  technique  assumes  the  ap  and 

rn 

aP  are  the  same  for  all  events  and  hence  that  each  event  is  at  the  same  depth. 

g 

The  slopes  of  the  calculated  travel-time  curves  imply  that  Vp^  =  6.021  km/sec 

and  Vp  =  7.736  km/sec,  values  which  are  to  be  compared  with  those  given  by 
n 

the  least-squares  fit  (Appendix  I)  of  Vp  =  6.036  km/sec  and  Vp  «  7.745  km/sec. 

g  n 
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TABLE  XX 


Absolute  errors  (km)  resulting  from  the  application  of  two  location  techniques 

to  a  common  data  set. 


Successive 

Determinations 

Pn  only 

Pn  and  Pg 

Pn  only 

P_  and  P 
n  g 

Event  II 

Pn 

II  P 

g 

Depth  Free 

Depth  Free 

Depth  Restrained 

Depth  Restrained 

FAULTLESS 

5 

5 

18.64 

(3.95) 

18.64 

4.05 

RULISON 

10 

7 

restr . 

restr. 

2.65 

3.38 

PASSAIC 

4 

3 

(2.66) 

restr. 

2.66 

2.04 

2.14 

a  3.75  * 

ROCKVILLE  DAM 

9 

8 

restr. 

restr. 

b  5.00 

DORMOUSE' 

3 

4 

diverges 

(9.39) 

diverges 

10.70 

KLICKITAT 

13 

12 

diverges 

(3.11) 

4.34 

4.11 

BANDICOOT 

5 

4 

4.18 

restr. 

4.18 

4.37 

SHOAL 

12 

0 

diverges 

no  Pg 

2.34 

no  P 

R 

MERRIMAC 

8 

0 

restr . 

no  Pg 

5.17 

no  P 

R 

CASBUGCY 

12 

12 

10.51 

restr. 

10.51 

10.02 

PILEDRIVER 

5 

5 

diverges 

(1.83) 

10.41 

2.56 

ROANOKE 

4 

4 

(4.41) 

ignores  Pg 

4.38 

4.24 

SALMON 

5 

0 

diverges 

no  Pg 

33.00 

no  Pg 

GNOME 

18 

1 

restr. 

restr. 

7.10 

7.06 

Simultaneous  Inversion 

Event 

Pn  only 

pn  and  Pg 

J.E.D. 

FAULTLESS 

22.35 

18.46 

18.31 

RULISON 

1.78 

0.69 

PASSAIC 

2.66 

2.75 

1.45 

ROCKVILLE  DAM 

2.59 

4.42 

DORMOUSE’ 

insuf.  Pn 

8.25 

4.72 

KLICKITAT 

0.97 

0.63 

1.18 

BANDICOOT 

11.14 

8.67 

5.05 

SHOAL 

1.88 

no  Pg 

MERRIMAC 

5.17 

no  Pg 

5.31 

GAS BUGGY 

10.56 

10.07 

PILEDRIVER 

8.81 

5.67 

6.13 

ROANOKE 

4.00 

4.11 

7.52 

SALMON 

107.75 

no  Pg 

GNOME 

7.15 

7.15 

*  Solution  alternates  between  two  epicenters 


In  all  versions  of  the  program  LOCATION,  arrival  time  residuals  for  each 

seismic  phase  are  weighted  by  the  reciprocals  of  the  standard  deviations  which 

are  anticipated  for  residuals  of  that  phase.  We  have  heretofore  applied  weights 

of  Wp^  =  1/op  =  1.0  and  Wp  =  1/3.0.  Table  X  presents  the  results  of 

changing  these  weighting  factors  in  the  method  of  simultaneous  inversions.  In 

accordance  with  the  observed  scatter  of  data  points  about  the  least-squares  lines 

in  Figure  4o ,  we  have  tried  alternate  values  of  Wp^  —  1/op^  =  and 

Wp  =  1/1.58.  It  is  not  the  values  of  these  weights  but  rather  their  ratio 
g 

which  affects  the  location,  so  we  anticipate  that  the  locations  which  will  be 

subject  to  the  most  change  by  these  new  weighting  factors  will  be  those  which 

utilize  the  most  values  of  P  relative  to  the  number  of  P  values.  Table  X 

g  n 

shows  that  only  in  the  case  of  PILEDRIVER  does  the  change  in  weights  signifi¬ 
cantly  improve  the  location,  and  in  the  case  of  ROCKVILLE  DAM,  which  was  not 
one  of  the  eight  NTS  events  used  to  determine  the  new  standard  deviations,  the 
new  location  is  significantly  worse.  The  value  of  JED  appears  to  be  enhanced 
if  the  new  weighting  scheme  is  applied;  this  improvement  is  illusory,  however, 
since  in  fact  the  locations  failed  to  converge  in  this  case.  The  values  which 
are  shown  for  the  absolute  errors  are  thus  unstable,  and  hence  they  are 
sensitive  to  the  locations  chosen  as  a  first  approximation. 
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TABLE  X 


Absolute  errors  (km)  resulting  from  assigning  more  weight 

to  P  data  in  the  method  of  simultaneous  inversions 
n 


Op  =1.0 

rn 

o_  =  0.84 
“n 

Op  =1.0 

*n 

op  =  0.84 
*n 

Op  =  3.0 
rg 

op  =  1.58 
pg 

op  =  3.0 
pg 

op  =1.58 
8 

Event 

No.  P 

n 

No.  P 

g 

Pn  and  Pg 

Pn  and  P 
n  g 

J.E.D. 

J.E.D.2 

FAULTLESS 

5 

5 

18.46 

21. 051 

18.31 

13.35 

RULISON 

10 

7 

0.69 

0.83 

PASSAIC 

4 

3 

2.75 

2.84 

1.45 

2.61 

ROCKVILLE  DAM 

9 

8 

4.42 

8.63 

DORMOUSE' 

3 

4 

8.25 

9.21 

4.72 

7.41 

KLICKITAT 

13 

12 

0.63 

0.32 

1.18 

1.00 

BANDICOOT 

5 

4 

8.67 

6.05 

5.05 

4.25 

MERRIMAC 

8 

0 

5.31 

3.78 

GAS BUGGY 

12 

12 

10.07 

9.66 

PILEDRIVER 

5 

5 

5.67 

2.48 

6.13 

1.03 

ROANOKE 

4 

4 

4.11 

4.30 

7.52 

5.94 

GNOME 

18 

1 

7.15 

7.15 

2  Pg  data  points  deleted 

failed  to  converge;  values  represent  the  minimum  found  for  x2 


2 


4.3  Error  Ellipses 


Thus  far  we  have  discussed  only  the  absolute  error  resulting  from  applying 
the  method  of  simultaneous  inversions  to  fourteen  events  whose  hypocenters  and 
origin  times  were  known  (X  pfviofu..  We  now  address  the  question  which  would  be 
encountered  in  the  location  of  unknown  events,  namely  the  a  pObt&U-OfU,  esti¬ 
mation  of  the  location  error  using  only  the  observed  data.  In  the  case  of 
conventional  location  determination  used  fixed  travel-time  relations,  the 
procedure  for  this  error  estimation  is  well  known  (Flinn,  1965).  Abbreviating 
the  system  of  equations  upon  which  the  location  algorithm  is  based  as 


5  c  =  B  (x  -  xQ)  , 


(29) 


where  6  is  the  vector  of  arrival  time  residuals,  x  is  the  unknown  location,  and 

-»■  ^  th  -*• 

Xq  is  the  trial  location,  the  standard  error  of  the  i  component  of  x  is 

given  by 


s  =  [(bt  B)“V/2 

ii 


(30) 


and  the  "error  ellipsoid"  surrounding  the  computed  hypocenter  x  is  given  by 


-  x)T  B_1(x  -  x)  <  N 


a  2  p 

0  N,N-4:0.05 


(31) 


where  N  is  the  number  of  dimensions  of  x,  where  E  is  an  NxN  submatrix  of 
(R  B)  ,  and  where  is  the  weighted  variance  of  the  final  residuals. 

We  see  that  the  method  of  simultaneous  inversions  presents  no  new  dif¬ 
ficulties  for  the  error  estimation,  since  it  can  operate  in  eight  dimen¬ 
sions  as  well  as  in  four.  In  Table  XI  we  present  the  standard  errors 
corresponding  to  the  absolute  location  errors  already  presented.  On 
account  of  the  previously  discussed  indetermination  of  the  origin  time,  5 ap 

“g 

is  understood  to  mean  the  standard  error  of  (T  +  ap  ),  and  similarly  for 

5  apn*  The  standard  errors  in  the  slopes  of  the  travel-time  curves  bp  and 

bp  have  been  used  to  compute  the  errors  in  their  reciprocals,  the  vefocities 
“n 

Vp  and  Vp  .  By  setting  the  number  of  dimensions  N  equal  to  2,  we  may 
calculate  the  error  ellipses  surrounding  the  computed  epicenters;  Figures  8a-n 
show  these  ellipses  along  with  the  actual  epicenters. 

For  those  events  for  which  the  degrees  of  freedom  (=  number  of  arrival 

times  -  number  of  parameters  solved  for  by  inversion)  are  small,  the  drror 

ellipses  are  unrealistically  large.  This  situation  can  be  alleviated  by 

using  xZ>  rather  than  F,  statistics  in  the  computation  (Evernden,  1969). 
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Sx  *  STANDARD  ERROR  C EAST  3  ©  TRUE 

Sr  *  STANDARD  ERROR  C NORTH  3  I - 

ELLIPSE  IS  FOR  CHI-SQUARED  STATISTIC 


EPICENTER 
H  8. 25M  KM 


Figure  8a  -  Error  Elipses  for  FAULTLESS 


S*  «  STANDARD  ERROR  C  EAST  3  ©  TRUE 

St  .  STANDARD  ERROR  C NORTH  3  I - 

OUTER  ELLIPSE  .  F-STATISTIC 

INNER  ELLIPSE  .  CHI -SQUARED  STATISTIC 


EPICENTER 
H  5* 729  KM 


Figure  8b  -  Error  Ellipses  for  RULISON 


2.  W*S* 


Sx  I  STANDARD  ERROR  CEASTJ  ©  TRUE  EPICENTER 

Sr  i  STANDARD  ERROR  C NORTH)  I - 1  7.713  KM 

ELLIPSE  IS  FOR  CHI-SQUARED  STATISTIC 


Figure  8c  -  Error  Ellipses  for  PASSAIC 
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Sx  .  STANDARD  ERROR  C EAST 3  ffi  TRUE  EPICENTER 

Sr  .  STANDARD  ERROR  C  NORTH  3  I - 1  5*603  KM 

OUTER  ELLIPSE  .  F-STATISTIC 

INNER  ELLIPSE  *  CHI-SQUARED  STATISTIC 


Figure  8d  -  Error  Ellipses  for  ROCKVILLE  DAM 


S*  «  STANDARD  ERROR  (EAST)  ©  TRUE  EPICENTER 

Sr  «  STANDARD  ERROR  (NORTH 3  I - 1  69*665  KM 

OUTER  ELLIPSE  »  F-STATISTIC 

INNER  ELLIPSE  »  CHI-SQUARED  STATISTIC 

Figure  8e  -  Error  Ellipses  for  DORMOUSE' 
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9* 757#S 


2.  480*Sr 


S»  .  STANDARD  ERROR  t  EAST  D  ®  TRUE  EPICENTER 

Sr  .  STANDARD  ERROR  C  NORTH  I  I - 1  4.299  KM 

OUTER  ELLIPSE  .  F-STATISTIC 

INNER  ELLIPSE  <  CHI-SQUARED  STATISTIC 

Figure  8f  -  Error  Ellipses  for  KLICKITAT 


2.  W*S* 


S*  «  STANDARD  ERROR  (EAST)  ©  TRUE  EPICENTER 
Sr  .  STANDARD  ERROR  (NORTH)  I - 1  16*151  KM 

ELLIPSE  IS  FOR  CHI -SQUARED  STATISTIC 


Figure  8g  -  Error  Ellipses  for  BANDICOOT 


S  *lbb'Z 


St  «  STANDARD  ERROR  C  EAST  3  ®  TRUE 

St  •  STANDARD  ERROR  C NORTH  3  I - 

OUTER  ELLIPSE  .  F-STATISTIC 

INNER  ELLIPSE  .  CHI -SQUARED  STATISTIC 


EPICENTER 
H  6. 309  KM 


Figure  8h  -  Error  Ellipses  for  SHOAL 


2. HH7*Sx 


Sx  «  STANDARD  ERROR  (EAST 3  ffi  TRUE  EPICENTER 

St  .  STANDARD  ERROR  (NORTH 3  I - 1  9*659  KM 

ELLIPSE  IS  FOR  CHI -SQUARED  STATISTIC 


Figure  81  -  Error  Ellipses  for  MERRIMAC 


S *Lbb*Z 


2» HH7*S* 


Sx  •  STANDARD  ERROR  (EAST 3  ffi  TRUE  EPICENTER 

St  •  STANDARD  ERROR  (NORTH!  I - 1  3,9HH  KM 

ELLIPSE  IS  FOR  CHI -SQUARED  STATISTIC 


Figure  8j  -  Error  Ellipses  for  GASBUGGY 
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Si  «  STANDARD  ERROR  (EAST 3  ffi  TRUE 

Sr  i  STANDARD  ERROR  (NORTH)  I - 

ELLIPSE  IS  FOR  CHI-SQUARED  STATISTIC 


EPICENTER 
H  8*  92*4  KM 


Figure  8k  -  Error  Ellipses  for  PILEDRIVER 


2*  W*S* 


8*  556*Sx 


Sx  .  STANDARD  ERROR  (EAST)  ©  TRUE  EPICENTER 

Sr  i  STANDARD  ERROR  (NORTH)  I - 1  305*883  KM 

OUTER  ELLIPSE  *  F-STATISTIC 

INNER  ELLIPSE  i  CHI-SQUARED  STATISTIC 


Figure  8m  -  Error  Ellipses  for  SALMON 
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2.  W*Si 


S«  .  STANDARD  ERROR  t EAST )  ffi  TRUE  EPICENTER 

Sr  •  STANDARD  ERROR  (NORTH)  I - 1  3.9t9  KM 

ELLIPSE  IS  FOR  CHI -SQUARED  STATISTIC 

Figure  8n  -  Error  Ellipses  for  GNOME 
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TABLE  XI 


Standard  Errors  -  Method  of  Simultaneous  Inversions 


Event 

(km) 

(km) 

“  r  n 
(sec) 

(sec/3eg) 

(secf 

(sec/aeg) (km/sec) 

(km/sic 

FAULTLESS 

5.36 

8.43 

1.0 

0.198 

2.8 

0.496 

0.104 

0.152 

RULISON 

4.14 

3.42 

1.1 

0.148 

3.2 

0.457 

0.089 

0.154 

PASSAIC 

7.69 

7.88 

1.6 

0.599 

5.8 

1.864 

0.337 

0.737 

ROCKVILLE 

DAM 

5.52 

4.54 

2.0 

0.302 

7.4 

1.074 

0.182 

0.326 

DORMOUSE ' 

8.82 

17.85 

3.0 

0.922 

4.3 

1.519 

0.530 

0.518 

KLICKITAT 

4.33 

2.89 

0.7 

0.144 

2.1 

0.391 

0.078 

0.130 

BANDICOOT 

4.39 

16.50 

3.5 

1.161 

6.8 

2.125 

0.746 

0.684 

SHOAL 

4.31 

4.04 

0.8 

0.142 

0.078 

MERRIMAC 

4.41 

4.76 

1.1 

0.248 

0.133 

GASBUGGY 

3.50 

4.03 

1.2 

0.167 

3.2 

0.441 

0.101 

0.150 

PILEDRIVER 

4.15 

8.49 

2.5 

0.785 

6.8 

2.180 

0.400 

0.663 

SALMON 

26.54 

89.38 

8.8 

0.223 

0.138 

GNOME 

2.52 

4.03 

0.7 

0.133 

0.076 
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If  we  assume  that  the  a  pfii.0fu.  assignation  of  standard  deviations  for  P 

and  P  arrival  time  measurements  are  correct  (  which,  as  we  have  shown,  they 
n 

are  not),  then  equation  (31)  should  be  replaced  by: 

<*  "  X')T  if1  (x  “  x)  -x22:0.05*  (32) 

In  this  equation  we  have  assumed  that  the  weighted  standard  deviations, 

°i^w±'  are  e<lua^  to  unity.  On  each  error  ellipse  we  show  the  effect  of 
changing  the  length  of  the  axes  in  this  way.  If  the  number  of  degrees  of 
freedom  of  the  solution  is  small,  then  may  be  an  underestimate  of  the  true 
variance  of  o^,  and  the  right-hand  side  of  equation  (31)  would  be  spuriously 
low.  In  those  cases  in  which  such  an  underestimate  results  in  a  smaller 
value  for  the  right-hand  side  of  equation  (31)  than  is  found  for  the  right- 
hand  side  of  equation  (32),  we  have  ignored  the  anomalously  small  F-statistic 
ellipse  and  have  plotted  only  the  ellipse  for  the  chi- squared  statistic. 


cUNCLUblONb  AND  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 


0  . 

Significant  improvements  have  been  made  in  locating  seismic  event 
hypocenters  with  the  successive  determinations  method  and  also  simultaneous 
inversions  method.  We  have  learned  that  the  earth  is  typically  hetero¬ 
geneous  and  that  no  simple  'model'  of  an  area  is  adequate  for  use  in  each 
single  case.  Location  experiments  of  Chang  and  Racine  (1979)  showed  no 
improvements  in  location  errors,  probably  because  the  regional  models  they 
used  were  not  adequate.  Experiments  with  local  models  resulted  in  small 
improvements,  but  with  no  significant  difference  from  using  one  simple 
(Herrin)  model.  This  is  because  local  models  adequate  for  the  source  area 
may  not  be  adequate  for  the  receiver  area,  and  vice  versa. 

We  found  that  the  HYLO  method,  which  uses  average  velocity  through 
the  source  to  station  path,  made  no  significant  improvements,  possibly  because 
the  earth  is  so  heterogeneous  that  the  contour  map  of  velocity  which  we  have 
used  in  the  experiments  is  not  suitable  for  each  case. 

The  best  result  can  be  obtained  by  modifying  P  and  P^  velocities  for 
each  event,  whether  successively  or  simultaneously.  This  means  that  even 
for  two  events  a  few  tens  of  kilometers  apart,  Pg  and  P^  velocities  may  be 
significantly  different,  perhaps  due  to  regional  and  local  heterogeneities  and 
their  influence  on  the  characters  of  the  radiated  waves. 

Note  that  in  these  two  methods  we  are  not  trying  to  obtain  local  crustal 

structures  by  inversions.  The  residual  data  from  the  successive  and  simultaneous 

methods  may  however  contain  information  regarding  the  dip  of  the  Moho.  Thus,  it 

may  be  possible  to  extract  crustal  models  using  these  methods.  In  addition,  the 

successive  and  simultaneous  methods  should  be  extended  to  include  more 

regional  phases,  including  L^,  which  is  generally  the  highest  amplitude  regional 

phase  on  the  record.  The  two  methods  should  be  merged  with  the  program  which 

utilizes  teleseismic  data,  and  used  to  evaluate  location  accuracies  for  various 

combinations  of  regional  and  teleseismic  data.  Tests  should  be  carried  out 

on  data  other  than  from  the  WUS.  The  methods  should  be  tested  on  AI  data. 

Since  the  methods  will  utilize  the  phases  P  and  L  ,  whose  start  times  are 

g  g 

often  unclear,  an  investigation  of  the  precision  with  which  analysts  and 
automatic  detectors  can  pick  the  onset  of  the  various  phases  should  be 
evaluated.  Finally,  the  program  should  be  incorporated  into  the  Regional 
Event  Location  System  (RELS) . 
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APPENDIX  I 


Least-Squares  Fit  to  Figure  4 
a:  intercept  of  least-squares  line 

b:  slope  of  least-squares  line  (reciprocal  of  velocity) 
JED:  8  NTS  events  used  together  in  Figure  4o. 


FAULTLESS  P 

n 

TT-Linear 

Fit 

Station 

A* 

TT  (sec) 

resid  (sec) 

JED  resid 

MN-NV 

1.528 

28.4 

0.14 

0.62 

KN-UT 

3.128 

50.6 

-0.71 

-0.15 

UBO 

5.412 

86.1 

1.89 

2.56 

TFO 

5.905 

90.1 

-1.21 

-0.51 

LC-NM 

ap  =  6.249  sec 

10.000 

150.2 

-0.11 

0.80 

bp^  =  14.406  sec/deg 

V  = 

7.719  km/sec 

FAULTLESS  P 

g 

Station 

A° 

TT  (sec) 

resid  (sec) 

JED  resid 

MN-NV 

1.528 

29.1 

0.31 

0.50 

KN-UT 

3.128 

58.2 

-0.83 

0.13 

UBO 

5.412 

102.1 

-0.11 

1.95 

TFO 

5.905 

112.4 

0.87 

3.17 

LC-NM 

10.000 

188.7 

-0.24 

4.04 

ap  =  -0.100  sec 

bpn  *  18.904  sec/deg 
*8 

V  = 

5.882  km/sec 

RULISON  P 

n 

Station 

A° 

TT  (sec) 

resid 

UBO 

1.546 

29.4 

-0.86 

KN-UT 

4.519 

71.2 

0.59 

ALQ 

4.612 

72.9 

1.03 

TFO 

5.778 

88.2 

0.50 

LC-NM 

7.074 

106.4 

1.11 

LAO 

7.388 

108.6 

-0.95 

TUC 

7.445 

108.6 

-1.72 

CR2NB 

8.607 

125.4 

-0.70 

BP-CL 

8.679 

125.9 

-1.17 

BMO 

ap  *  9.279  sec 

8.823 

131.2 

2.17 

bpn  *  13.572  sec/deg  ■+• 

V  = 

8.192  km/sec 
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RULISON  P 

g 

Station 

A° 

APPENDIX  1 

(Continued) 

XT 

resid 

UBO 

1.546 

31.3 

2.22 

KN-UT 

4.519 

80.7 

-2.47 

TFO 

5.778 

103.4 

-2.68 

LC-NM 

7.074 

130.9 

1.24 

CR2NB 

8.607 

159.1 

1.55 

BP-CL 

8.679 

159.9 

1.04 

BMO 

8.823 

160.6 

-0.89 

ap  =  0.950  sec 

bp|  =  18.195  sec/deg  V  =  6.11  km/sec 


PASSAIC  P 

n 


Station 

A° 

TT 

resid 

JED  resid 

DV-CL 

1.274 

24.5 

0.04 

0.37 

MN-NV 

2.125 

36.3 

-0.27 

-0.05 

KN-UT 

2.571 

43.2 

0.28 

0.45 

FS-AZ 

4.343 

68.1 

-0.04 

-0.09 

» 

ap  =  6.330 

sec 

bp^  =  14.232  sec/deg  V  = 

7.812  km/sec 

% 

PASSAIC  P 

g 

Station 

A° 

TT 

resid 

JED  resid 

MN-NV 

2.125 

39.7 

-0.36 

0.10 

I 

KN-UT 

2.571 

48.1 

0.45 

0.29 

1 

FS-AZ 

4.343 

77.7 

-0.09 

-2.76 

ap  =  3.913  sec 

bpg  =  17.011  sec/deg  -+■  V  =  6.536  km/sec 
g 


ROCKVILLE  DAM  P 

n 


Station  A0 

TT 

resid 

A2 

TT2 

UBO  2.579 

46.0 

0.41 

6.651 

2116.0 

KN-UT  5.531 

85.7 

0.14 

30.591 

7344.5 

WN-SD  6.112 

92.7 

-0.72 

37.356 

8593.3 

RG-SD  6.241 

95.2 

0.03 

38.950 

9063.0 

TFO  6.380 

97.5 

0.45 

40.704 

9506.3 

LAO  7.327 

109.2 

-0.67 

53.685 

11924.6 

CR-NB  7.488 

111.9 

-0.15 

56.070 

12521.6 

WMO  7.814 

115.8 

-0.67 

61.058 

13409.6 

KC-MO  9.133 

135.5 

1.18 

83.412 

18360.3 

ap  "  10.676  sec 

bpn  ■  13.538  sec/deg 
rn 

->•  V  =  8. 

213  km/sec 
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ROCKVILLE  DAM  P 

g 

Station  A0  TT  resid 


KN-UT 

5.531 

101.3 

-1.81 

WN-SD 

6.112 

114.2 

0.14 

RG-SD 

6.241 

116.9 

0.41 

TFO 

6.380 

117.1 

-2.01 

LAO 

7.327 

143.5 

6.55 

CR-NB 

7.488 

140.0 

0.01 

WMO 

7.814 

145.0 

-1.13 

KC-MO 

9.133 

168.8 

-2.19 

aP„  = 

-1.129 

sec 

bp®  = 

18.846 

sec/deg  V  = 

5.900  km/sec 

DORMOUSE  '  P 

n 


Station 

A° 

TT 

resid 

JED  resid 

MN-NV 

2.188 

37.4 

-0.13 

0.15 

KN-UT 

2.555 

42.9 

0.16 

0.38 

FS-AZ 

4.293 

67.4 

-0.03 

-0.07 

aPn  = 

6.454 

sec 

bPn« 

14.203 

sec/deg 

-*  V  = 

7.829  km/ sec 

KLICKITAT  P 

n 

Station 

A° 

TT 

resid 

JED  resid 

CU-NV 

1.592 

28.6 

0.09 

-0.10 

EK-NV 

2.071 

35.9 

0.53 

0.33 

MN-NV 

2.106 

36.2 

0.32 

0.12 

KM- CL 

2.470 

41.3 

0.21 

0.00 

KN-UT 

2.572 

42.8 

0.25 

0.03 

CP-CL 

4.420 

68.2 

-0.82 

-1.10 

TFO 

4.831 

74.8 

-0.11 

-0.40 

BX-UT 

5.279 

80.2 

-1.13 

-1.43 

UBO 

5.968 

93.3 

2.10 

1.78 

DR-CO 

6.580 

97.9 

-2.07 

-2.40 

HL-ID 

6.631 

99.5 

-1.20 

-1.54 

PI-WY 

7.279 

110.3 

0.32 

-0.04 

LC-NM 

9.095 

137.5 

1.51 

1.09 

aPn  ’ 

5.708 

sec 

bPn  * 

14.325 

sec/deg 

+  V  « 

7. 762  km/sec 
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(Continued) 


KLICKITAT  P 

g 

Station  A° 

TT 

resid 

JED  resid 

CU-NV  1.592 

30.5 

0.80 

0.72 

EK-NV  2.071 

38.4 

-0.04 

-0.20 

KM-CL  2.470 

46.6 

0.88 

0.65 

KN-UT  2.572 

47.0 

-0.59 

-0.83 

CP-CL  4.420 

79.7 

-1.62 

-2.17 

TFO  4.831 

89.7 

0.88 

0.26 

BX-UT  5.279 

96.7 

-0.30 

-1.00 

UBO  5.968 

110.7 

1.12 

0.31 

DR-CO  6.580 

119.7 

-1.05 

-1.96 

HL-ID  6.631 

120.7 

-0.99 

-1.90 

PI-WY  7.279 

ap  =  0.635  sec 

132.9 

1.54 

0.21 

bpg  =  18.255  sec/deg  -*•  V  = 
g 

6.091  km/ sec 

BANDICOOT  P 

n 

Station  A0 

TT 

resid 

JED  resid 

MN-NV  2.188 

37.4 

-0.29 

0.15 

KN-UT  2.555 

43.0 

0.28 

0.48 

TF-CL  3.704 

58.7 

0.21 

-0.32 

FS-AZ  4.293 

66.2 

-0.38 

-1.27 

CP-CL  4.311 

ap  =  7.656  sec 

67.0 

0.18 

-0.73 

bPn  =  13.725  sec/deg  -+  V  = 
rn 

8.101  km/sec 

BANDICOOT  P 

g 

Station  A° 

TT 

resid 

JED  resid 

MN-UT  2.188 

40.5 

-0.04 

-0.26 

KN-UT  2.55 

48.0 

0.38 

0.48 

TF-CL  3.704 

68.8 

-0.99 

0.12 

FS-AZ  4.293 

81.8 

0.65 

2.27 

ap  =  -1.680  sec 

bPg  =  19.295  sec/deg  •+  V 

=  5.763  km/sec 

g 
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SHOAL  P 

n 

Station 

A° 

TT 

resid 

EK-NV 

2.074 

35.6 

0.25 

WI-NV 

2.261 

37.8 

-0.21 

MV-CL 

2.263 

37.7 

-0.34 

CU-NV 

2.340 

39.3 

0.17 

KN-UT 

4.889 

75.2 

-0.13 

HL-ID 

5.420 

81.5 

-1.37 

BMO 

5.701 

87.8 

0.95 

CP-CL 

6.659 

100.2 

-0.26 

UBO 

6.881 

106.3 

2.69 

BX-UT 

7.215 

107.5 

-0.85 

TFO 

7.536 

113.0 

0.09 

DR-CO 

8.498 

125.6 

-0.97 

3pn  *  5.906  sec 

bp  *  14.199  sec/deg  -*•  7.831  km/sec 
n 


MERRIMAC  P 

n 


Station 

A° 

TT 

resid 

JED  resid 

MN-NV 

2.170 

37.2 

0.56 

0.21 

KN-UT 

2.565 

42.6 

0.25 

-0.07 

FS-AZ 

4.309 

67.8 

0.23 

0.10 

CP-CL 

4.325 

66.8 

-1.01 

-1.13 

WI-NV 

4.430 

68.8 

-0.52 

-0.64 

MV-CL 

4.673 

72.6 

-0.24 

-0.33 

DR-CO 

6.584 

100.3 

-0.17 

-0.06 

HL-ID 

6.724 

103.4 

0.90 

1.03 

ap  =  5.261  sec 

bp^  =  14.461  sec/deg  ■>  7.689  km/sec 


GASBUGGY  P 

n 


Station 

A° 

TT 

resid 

TFO 

4.029 

65.6 

0.30 

UBO 

4.047 

66.1 

0.56 

LC-NM 

4.305 

67.4 

-1.61 

KN-UT 

4.440 

71.7 

0.87 

PQ-ID 

6.507 

98.4 

-0.26 

WMO 

7.359 

108.6 

-1.52 

WZ-NV 

7.402 

111.9 

1.20 

WN-SD 

8.537 

125.4 

-0.58 

HL2ID 

8.768 

130.7 

1.61 

MN-NV 

8.786 

129.7 

0.37 

BS-MA 

9.124 

134.6 

0.72 

GV-TX 

ap  *  11.064 
bP"  -  13.461 

9.285 

sec 

sec/deg  V  ■ 

134.4 

8.260  km/sec 
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GASBUGGY  P 

g 


Station 

A°  . 

TT 

resid 

TFO 

4.029 

75.2 

-0.77 

UBO 

4.047 

76.1 

-0.20 

LC-NM 

4.305 

79.1 

-1.85 

KN-UT 

4.440 

86.9 

3.51 

PQ-ID 

6.507 

121.2 

0.54 

WMO 

7.359 

132.6 

-3.42 

WZ-NV 

7.402 

138.4 

1.61 

WN-SD 

8.537 

157.4 

0.14 

HL2ID 

8.768 

159.6 

-1.83 

MN-NV 

8.786 

161.1 

-0.65 

BS-MA 

9.124 

170.9 

3.06 

GV-TX 

9.285 

170.6 

-0.15 

ap  =  3.324 

sec 

bPg  =  18.032 

sec/deg  -+  V  = 

6.166  km/ 

sec 

PILEDRIVER  P 

n 

Station 

A° 

TT 

resid 

JED  resid 

MN-NV 

2.049 

35.2 

-0.10 

-0.06 

KN-UT 

2.589 

43.1 

0.03 

0.09 

TFO 

4.885 

75.7 

-0.40 

-0.27 

UBO 

5.936 

92.1 

0.88 

1.04 

BMO 

7.673 

115.8 

-0.41 

-0.20 

ap  =  5.823 

sec 

bp"  =  14.386 
rn 

sec/deg  -*■  V  = 

7.729  km/sec 

PILEDRIVER  P 

g 

Station 

A0 

TT 

resid 

JED  resid 

MN-NV 

2.049 

37.6 

-0.14 

-0.60 

KN-UT 

2.589 

47.0 

-0.53 

-1.15 

TFO 

4.885 

90.6 

1.41 

0.16 

UBO 

5.936 

108.0 

-0.26 

-1.80 

BMO 

7.673 

139.3 

-0.48 

-2.49 

ap  =  0.558 

sec 

bpg  =  18.144 
g 

sec/deg  •+  V  = 

6.128  km/ 

sec 

ROANOKE  Pn 

Station 

A0 

TT 

resid 

JED  resid 

MN-NV 

2.116 

36.1 

0.04 

-0.12 

KN-UT 

2.580 

42.9 

-0.04 

0.02 

FM-UT 

3.681 

58.9 

-0.38 

0.12 

TF-CL 

3.724 

60.3 

0.38 

1.00 

ap  *=  4.662 

sec 

bpn  =  14.838 

sec/deg  -*■  V  = 

7.494  km/ 

sec 

rn 
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ROANOKE  Pg 


Station 

A° 

TT 

resid 

JED  resid 

MN-NV 

2.116 

40.7 

0.77 

1.27 

KN-UT 

2.580 

47.8 

-1.05 

-0.18 

FM-UT 

3.681 

69.1 

-0.90 

0.84 

TF-CL 

3.724 

72.0 

1.18 

2.95 

bp®  =  19.212  sec  deg  5.788  km/sec 
rg 

SALMON  P 

n 


Station 

A° 

TT 

resid 

EU-AL 

2.178 

36.5 

-0.29 

JE-LA 

2.185 

37.4 

0.52 

CPO 

5.560 

82.4 

0.11 

GV-TX 

6.544 

94.6 

-0.94 

WMO 

8.377 

120.8 

0.60 

aP  =  7.482 

sec 

bPn  =  13.456 
rn 

sec/deg  -*■  V 

=  8.263  km/sec 

GNOME  P 

n 

Station 

A° 

TT 

resid 

LC-NM 

2.319 

40.6 

0.70 

PO-TX 

2.429 

42.1 

0.65 

SS-TX 

2.594 

43.8 

0.03 

SM-TX 

4.177 

65.6 

-0.39 

ML-NM 

4.339 

68.2 

-0.06 

RT-NM 

4.474 

70.2 

0.04 

LP-TX 

4.738 

73.0 

-0.87 

SV-AZ 

4.820 

75.5 

0.48 

HB-OK 

5.049 

78.1 

-0.13 

AM-OK 

5.698 

85.7 

-1.64 

SF-AZ 

5.972 

91.5 

0.31 

DR-CO 

6.112 

94.6 

1.45 

TO-OK 

6.456 

95.6 

-2.38 

SJ-TX 

6.685 

99.7 

-1.50 

FS-AZ 

6.808 

103.7 

0.77 

MW-AZ 

7.625 

116.6 

2.20 

KN-UT 

8.777 

131.2 

0.63 

PM-WY 

9.007 

133.5 

-0.30 

ap  “  7.348 

sec 

bD  14.039  sec/deg  -*■  V  ■ 

7.920  km/sec 

GNOME  P 

g 

Station 

A° 

TT 

LC-NM 

2.319 

44.0 
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APPENDIX  I 
(Continued) 

JED:  FAULTLESS,  PASSAIC,  DORMOUSE  PRIME,  KLICKITAT,  BANDICOOT,  MERRIMAC, 

PILEDRIVER,  and  ROANOKE  combined 


P^  (47  data  points) 

ap  =  5.842  sec 
n 

b„  =  14.356  sec/deg  V  =  7.745  km/sec 

Ml 

Op  =  0.835  sec 
n 

Pg  (37  data  points) 

ap  =  0.455  sec 
g 

bp  =  18.420  sec/deg  -*■  V  =  6.036  km/sec 

g 

Op  =  1.575  sec 
g 


» 
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